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Abstract 

The normal/superconducting phase boundary Tc has been calculated for meso- 
scopic loops, as a function of an applied perpendicular magnetic field H. While for 
thin- wire loops and filled disks the Tc{H) curves are well known, the intermedi- 
ate case, namely mesoscopic loops of finite wire width, have been studied much 
less. The linearized first Ginzburg-Landau equation is solved with the proper nor- 
mal/vacuum boundary conditions both at the internal and at the external loop 
radius. For thin-wire loops the Tc{H) oscillations are perfectly periodic, and the 
Tc{H) background is parabolic (this is the usual Little-Parks effect). For loops of 
thicker wire width, there is a crossover magnetic field above which Tc{H) becomes 
quasi-linear, with the period identical to the Tc{H) of a filled disk (i.e. pseudoperi- 
odic oscillations). This dimensional transition is similar to the 2D-3D transition for 
thin films in a parallel field, where vortices start penetrating the material as soon as 
the film thickness exceeds the temperature dependent coherence length by a factor 
1.8. For the presently studied loops, the crossover point is controlled by a similar 
condition. In the high field '3D' regime, a giant vortex state establishes, where only 
a surface superconducting sheath near the sample's outer radius is present. 
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1 Introduction 



The nucleation of superconductivity in mesoscopic samples has received a re- 
newed interest after the development of nanofabrication techniques, hke elec- 
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tron beam lithography. A superconductor is in the mesoscopic regime when 
the sample size is comparable to the superconducting coherence length ^(T). 
In the framework of the Ginzburg-Landau (GL) theory, ^(T) sets the length 
scale for spatial variations of the modulus of the superconducting order param- 
eter I \1' I . The pioneering work on mesoscopic superconductors was carried out 
already in 1962 by Little and Parks [1], who measured the shift of the critical 
temperature Tc{H) of a (multiply connected) thin- walled Sn microcylinder (a 
thin-wire "loop") in an axial magnetic field H. The Tc{H) phase boundary 
showed a periodic component, with the magnetic period corresponding to the 
penetration of a superconducting flux quantum $0 = h/2e. 

A few years later. Saint- James calculated the Tc{H) of a singly-connected 
cylinder [2] (a mesoscopic "disk"). Taking into account the analogy with the 
situation of a semi-infinite superconducting slab in contact with vacuum [3], 
the critical field was called H^ziT) in this case, since superconductivity nucle- 
ates initially near the sample interface. In the present paper, we will use the 
notation H*^{T), for the nucleation magnetic field. 

The Tf.{H) phase boundary (or H*^{T)) of the disk shows, just like for the 
usual Little-Parks effect in a multiply connected sample (loop), an oscillatory 
behaviour. When moving along Tc{H), superconductivity concentrates more 
and more near the sample interface as H grows. A giant vortex state is formed: 
a "normal" core carries L flux quanta, and the 'effective' loop radius increases, 
resulting in a decrease of the magnetic oscillation period. An experimental 
verification of these predictions was carried out later on by Buisson et al. [4] 
and by Moshchalkov et al. [5]. 

In the early paper from Saint- James and de Gennes [3], H*^{T) has been 
calculated also for a film exposed to a parallel magnetic field, where surface 
superconductivity can grow along two superconductor /vacuum interfaces. For 
low magnetic fields, the two surface superconducting sheaths overlap, and, 
as a result, Tc versus H becomes parabolic, which is characteristic for 2D 
behaviour. When increasing the field, a crossover to a linear Tc{H) dependence 
(3D) occurs at t ~ 2 ^(T), with t the film thickness. Shortly after, it was shown 
that vortices start to nucleate in the film at this dimensional crossover point 
(t = 1.8e(T))[6]. 

The goal of the present paper is to study the phase boundary Tc{H) of 
loops made of finite width wires. In a Type-II material, superconductivity 
is expected to be enhanced, with respect to the bulk upper critical field 
Hc2 '. H*^{T) > Hc2{T), both at the external and the internal sample sur- 
faces. As for a film in a parallel field, a 2D-3D dimensional crossover can 
be anticipated, since the loops may be simply considered as a film, which is 
bent such that its ends are joined together. We calculate, for the first time, 
the phase boundary Tc{H) as the ground state solution of the linearized first 



GL equation with two superconductor /vacuum interfaces. 

2 The linearized GL equation 

The hnearized first GL equation to be solved in order to find Tc{H) is: 

1 



2m* 



[-ihV - 2eA)^^ = I - a| ^ . 



This equation is formally identical to the Schrodinger equation for a particle 
with a charge 2e in a magnetic field. At the onset of superconductivity, the 
nonlinear GL term can be omitted and the z-dependence disappears from the 
equations and therefore an infinitely long cylinder and a disk have identical 
Tc{H) boundaries. It is further assumed that fioH = rotA, with H the applied 
magnetic field. The eigenenergies \ — a\ can be written as: 

/i^ _ ;i^ / T \ 
' ~ "' ~ 2m* e{T) ~ 2m* ^^(0) V 7^/ ' ^^ 



Tco being the critical temperature in zero magnetic field. From the energy 
eigenvalues of Eq. (1), the lowest Landau level | —aLLL{,H)\ is directly related 
to the highest possible temperature Tc{H), for which superconductivity can 
exist. 

For the loop geometries, we choose the cylindrical coordinate system (r, (f) 
and the gauge A = {noHr/2) e^, where e^ is the tangential unit vector. The 
exact solution of the Hamiltonian (Eq. (1)) in cylindrical coordinates takes 
the following form [7,8]: 



M/(<l>, V.) = e-^^ (^) exp [-^) Ki-n, L + 1, ^/%) (3) 

K{a, c, y) = ci M{a, c, y) + C2 U{a, c, y) . 

Here $ = fioHirr'^ is the applied magnetic flux through a circle of radius r. 
The number n determines the energy eigenvalue. Most generally, the function 
K{a, c, y) can be any linear combination of the two confluent hypergeometric 
functions (or Kummer functions) M{a,c,y) and U{a,c,y) [9], but the sam- 
ple topology puts a constraint on ci, C2, and n, via the Neumann boundary 
condition: 

{-ihV - 2eA)'^ = , (4) 



which the solutions ^ of Eq. (1) have to fulfill at the sample interfaces b. 
The eigenenergies of Eq. (1) can be written in the form: 

r„2 ri /, T,(H)\ , / 1\ * ,„„ * 



where $ = ^qHtttI is arbitrarily defined. The integer number L is the phase 
winding, or fluxoid quantum number. The parameter n depends on L and is 
not necessarily an integer number^ as we shall see further on. 

The bulk Landau levels are obtained when substituting n = 0,1,2, ■■■ in 
Eq. (5), meaning that the lowest level n = corresponds to the bulk upper 
critical field fioHc2iT) = %/ {2Tr^\T)). 

For a disk geometry [2,4], we have to take C2 = in Eq. (3) in order to avoid the 
divergency of U{a, c, y -^ 0) = oo at the origin. Selecting the lowest Landau 
level at each value $, one ends up with a cusp-like Tc{H) phase boundary [2], 
which is composed of values n < in Eq. (5), thus leading to H*^(T) > Hc2(T). 
A similar calculation was performed for a single circular microhole in a plane 
film ("antidot") [7], where ci = in Eq. (3), since M{a,c,y -^ oo) = oo. Here 
as well, the lowest Landau level consists of solutions with n < 0. At each cusp 
in Tc(<l>), the system makes a transition L ^ L±l, i.e. a fiux quantum enters 
or is removed from the sample. 

The loops we are currently studying have two superconducting/vacuum inter- 
faces, one at the outer radius To, and one at the inner radius r^. Consequently, 
the boundary condition (Eq. (4)) has to be fulfilled at both ro, and r^. As a 
result, we have a system of two equations and two variables n and C2 (ci = 1 
is chosen), which we solved for different values of ri/ro. 



3 Results 



Fig. 1 shows the Landau level scheme (dashed lines) calculated from Eqs. (3)- 
(5), for a loop with ri/ro = 0.5. The applied magnetic fiux $ = ^oHirrl is 
defined with respect to the outer sample area. The Tc{H) boundary is com- 
posed of \1' solutions with a different phase winding number L and is drawn 
as a solid cusp-like line in Fig. 1. At $ ~ 0, the state with L = is formed 
at Tc(<l>) and one by one, consecutive flux quanta L enter the loop as the 
magnetic field increases. For low magnetic fiux, the background depression of 
Tc is parabolic, whereas at higher fiux, Tc{^) becomes quasi-linear, just like 
for the case of a filled disk. The crossover point from parabolic to quasi-linear 
appears at about $ ~ 14 $o- 



Fig. 1. Calculated energy level scheme (dashed lines) for a superconducting loop 
with the ratio of inner to outer radius ri/ro = 0.5. The solid and dotted lines 
correspond to Hc2(T) and Hc3{T), respectively. 

Fig. 2. Inverse enhancement factor r/"^ = e{H*^ / e{H c2) (see Eq. (5)) for loops with 
different aspect ratio, compared to the case of a disk and an antidot. The horizontal 
dashed line at 77"^ = 0.59 = 1/1.69 corresponds to HcsiT) / Hc2(T) = 1.69 for a 
plane superconductor /vacuum boundary [3]. 

The solid and dotted straight lines in Fig. 1 are the bulk upper critical 
field Hc2{T) and the surface critical field Hc3{T) for a semi-infinite slab, 
respectively. In these units the slopes of the curves (see Eq. (5)) are e = 2 
for Hc2 (substitute n = in Eq. (5)) and e = 2/1.69 for Hcs- The ra- 
tio 1] = e{H C2) / ^{H cz) = 1-69 corresponds then to the enhancement factor 
Hc3(T)/Hc2(T) at a constant temperature. For the loops we are studying here, 
1] = t{H (.2) / ^{H*.^) is varying as a function of the magnetic field. 

The energy levels below the Hc2 line (solid straight line in Fig. 1) could be 
found by fixing a certain L, and finding the real numbers n < numeri- 
cally after inserting the general solution (Eq. (3)) into the boundary con- 
dition (Eq. (4)). Note that the lowest Landau level always has a lower en- 
ergy I — «($)! than for a semi-infinite superconducting slab, which implies 
Hl^{T)>H,3{T) = l.mH,2{T). 

As mentioned earlier, in a thin film of thickness t in a parallel field H, a 
dimensional crossover is found at t = 1.84 ^(T). For low fields (high ^) Tc{H) 
is parabolic (2D), and for higher fields vortices start penetrating the film 
and consequently Tc{H) becomes linear (3D) [6]. In Fig. 1 the small arrow 
indicates the point on the phase diagram Tc{^) where w = 1.84 ^(T). For the 
loops as well, the dimensional transition shows up approximately at this point, 
although the vortices are not penetrating the sample area in the 3D regime. 
Instead, the middle loop opening contains a coreless 'giant vortex' with an 
integer number of fiux quanta L <l>o- 

In order to compare the fiux periodicity of Tc(<l>), we have plotted, in Fig. 2, 
the lowest energy levels of Fig. 1 as r]"^ = e{H*^)/e{Hc2), for loops with a 
different ri/vo- In this representation, the dotted horizontal line at rj^^ = 0.59 
corresponds to the surface critical field line Hc3{T). The nucleation field of 
a disk H*^{T) > 1.69 Hc2{T), and for a circular microhole in an infinite film 
("antidot") [7] H*^{T) < 1.69H^2iT). As $ grows (the radius goes to infinity) 
the H*^(T) of both the disk and the antidot approaches the Hc3{T) line. 

For all the loops we study here, the presence of the outer sample interface 
automatically implies that H*^{T) > Hc3{T) is enhanced (77 > 1.69), with 
respect to the case of a fiat superconductor- vacuum interface. For loops with a 
small ri/vo, the Tc{^) boundary very rapidly collapses with the Tc{^) of the dot 



[r] becomes the same). The presence of the opening in the sample is not relevant 
for the giant vortex formation in the high flux 3D regime. On the contrary, in 
the low flux regime, the surface sheaths along the two interfaces overlap, giving 
rise to a different periodicity of Tc($) and to a parabolic background. This 
regime can be described within the London limit [10], since superconductivity 
nucleates almost uniformly within the sample. 

In summary, we have solved the linearized GL equation for loops of dif- 
ferent wire width, with Neumann boundary conditions at both the outer 
and the inner loop radius. The critical fields H*^{T) are always above the 
HcsiT) = 1.69 Hc2(T). The ratio H*^(T)/Hc2iT) is enhanced most strongly 
when the sample's surface-to-volume ratio is the largest. The Tc(<l>) behaviour 
can be split in two regimes: for low flux, the background of Tc is parabolic 
(characteristic for 2D behaviour) and the Little-Parks Tc{^) oscillations are 
perfectly periodic. In the high flux regime, the period of the Tc{^) oscillations 
is decreasing with $ and the background T^ reduction is quasi-linear (3D 
regime). The 2D-3D crossover between the two regimes, at a certain applied 
flux $, is similar to the dimensional transition in thin films subjected to a par- 
allel field. As soon the 3D regime is reached, a giant vortex state is created, 
where only a sheath close to the sample's outer interface is superconducting. 
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